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Abstract 



i-C . A method for quantum corrections of Hanbury-Brown/Twiss (HBT) interferometric radii pro- 

o ' duced by semi-classical event generators is proposed. These corrections account for the basic 

D ■ 

,^ ' indistinguishability and mutual coherence of closely located emitters caused by the uncertainty 



principle. A detailed analysis is presented for pion interferometry in p + p collisions at LHC energy 
(-y/s = 7 TeV). A prediction is also presented of pion interferometric radii for p + Pb collisions at 



QQ ' t/s = 5.02 TeV. The hydrodynamic/hydrokinetic model with UrQMD cascade as 'afterburner' is 

utilized for this aim. It is found that quantum corrections to the interferometry radii improve sig- 



nificantly the event generator results which typically overestimate the experimental radii of small 
systems. A successful description of the interferometry structure of p + p collisions within the 



k> ■ corrected hydrodynamic model requires the study of the problem of thermalisation mechanism, 

$-( < 

C^ . still a fundamental issue for ultrarelativistic ^4 + A collisions, also for high multiplicity p + p and 

p + Pb events. 

PACS numbers: 13.85.Hd, 25.75.Gz 



PACS: 24.lO.Nz, 24. 10. Pa, 25.75.-q, 25.75.Gz, 25.75.Ld. 

Keywords: correlation femtoscopy, HBT radii, proton-proton collisions, proton-nucleus colli- 
sions LHC, uncertainty principle, coherence 

I. INTRODUCTION 

The quantum-statistical enhancement of the pairs of identical pions produced with close 
momenta was observed first in p -\- p collisions in 1959 [1]. It took more than a decade 
to develop the method of pion interferometry based on the discovered phenomenon. This 
was done at the beginning of the 1970s by Kopylov and Podgoretsky 2|. Their theoretical 
analysis assumed the radiating source as consisting of independent incoherent emitters. In 
fact, such a representation is used for a long time for the analysis of the space-time structure 
of particle sources creating in p -\- p, p -\- p, e+ + e~ and A -|- A collisions. The concept of 



independent emitters was applied to a further development of the interferometric met 
in particular, to account for momentum-position correlations of the emitted particles 



lod, 

6| 



that, in turn, has resulted in a general interpretation of the measured radii as the homogene- 
ity lengths in the Wigner functions [7H9|. This concept is important for a study oi A -\- A 
collision processes within the hydrodynamic approach. Also a detailed analysis of the par- 
ticle final state (Coulomb) interactions brings the significant contribution to the traditional 



method of correlation femtoscopy 10|, lll |. 



In a recent paper 12|] the correlation analysis is taken beyond the model of independent 
particle emitters. It is found that the uncertainty principle leads to (partial) indistinguisha- 
bility of closely located emitters that fundamentally impedes their full independence and 
incoherence. The partial coherence of emitted particles is because of the quantum nature of 
particle emission and happens even if there is no specific mechanism to produce a coherent 
component of the source radiation. This effect leads to a reduction of the interferometry 
radii and suppression of the Bose-Einstein correlation functions. The effect is significant 



only for small sources with typical sizes less than 2 fm. We shall apply this approach 12 1 
to the analysis of data in p -\- p collisions at the LHC energy of ^/s = 7 TeV, where the 
measured interferometry radii are just within the above scale. A simple estimate will be 
done also for p + Pb, where the radii are larger and such corrections are less important. 
A first attempt of the systematic theoretical analysis of the pion interferometry oi p -\- p 



collisions at the top RHIC and a/s = 0.9 TeV LHC energies was made in Ref.[l3| within 
the Quark-Gluon String Model (QGSM). It was found that, for a satisfactory description of 
the interferometry radii, one needs to reduce significantly the formation time by increasing 
the string tension value relative to the one fixed by the QGSM description of the spectra 
and multiplicity. Otherwise, the radii obtained within QGSM are too large compared to 
the measured ones. Hypothetically one can hope to reduce the predicted radii suggesting 
the other approach - the hydrodynamic mechanism of the bulk matter production in p + p 
collisions, at least, for high multiplicity events. Then, to reproduce high multiplicity, the 
initially very small p+p system has to be superdense at early times. This leads to very large 
collective velocity gradients, and so the homogeneity lengths should be fairly small. How- 
ever, as we shall demonstrate, even at the maximally possible velocity gradients at the given 
multiplicity, one gets again an overestimate of the interferometry radii in p + p collisions. 
Therefore, one can conclude that the problem of theoretical description of the interferom- 
etry radii in p + p collisions may be a general one for different types of event generators 
associated with various particle production mechanisms. Here we try to correct the results 
on interferometry from event generators using for this aim the quantum effects accounting 
for partial indistinguishability and mutual coherence of the closely located emitters due to 
the uncertainty principle |12| . 



In this article we employ the hydrokinetic model (HKM) [1J| in its hybrid form [15 1 
where the UrQMD hadronic cascade is considered as the semi-classical event generator 
at the post freeze-out ("afterburner") stage of the hydrodynamic/hydrokinetic evolution. 
We analyze two aspects of the analysis oi p + p collisions. The main one is: whether 
quantum corrections can help to describe the experimental data. If yes, it gives hope that 
it can be successfully applied for any event generator associated with another mechanisms 
of the particle production. The second aspect is more sophisticated: whether the typical 
hybrid models developed for A + A collisions (here hybrid = hydrodynamic/hydrokinetic 
+ hadronic cascade) with correspondingly modified initial conditions and with the above- 
mentioned quantum corrections can be real candidates to describe the bulk observables in 
p + p collisions at LHC energies. For this aim we study the space-time structure oi p + p 
collisions, namely, analyze the multiplicity dependence of interferometry radii and volume as 
well as the p^-behavior of the HBT radii. It is worth noting that a satisfactory description 
of the corresponding experimental data challenges the theoretical picture oi p + p collisions. 



however supporting the Landau's pioneer suggestion [l6| to use relativistic hydro dynamic 
theory for the hadron colhsions with high multiphcity. 

II. HYDROKINETIC MODEL: DESCRIPTION AND RESULTS FOR p + p COL- 
LISIONS 

The hydrokinetic model [ij, [l5| oi A + A colhsions consists of several ingredients describ- 
ing different stages of the evolution of matter in such processes. At the first stage of system's 
evolution the matter is supposed to be chemically and thermally equilibrated and its expan- 



sion is described within perfect (2+l)D boost-invariant relativistic 



lydrodynamics with the 



lattice QCD-inspired equation of state in the quark-gluon phase [17| matched with a chemi- 
cally equilibrated hadron-resonance gas via crossover-type transition. The hadron-resonance 
gas consists of 329 well-established hadron state^j made of u,d,s-quarks, including a-meson 
(/o(600)). With such an equilibrated evolution the system reaches the chemical freeze-out 
isotherm with the temperature T^h = 165 MeV. At the second stage with T < Teh, the 
hydrodynamically expanding hadron system gradually looses its (local) thermal and chemi- 
cal equilibrium and particles continuously escape from the system. This stage is described 
within the hydrokinetic approach [ij] to the problem of dynamical decoupling. In hHKM 



model [15| the hydrokinetic stage is matching with hadron cascade UrQMD one [19| at the 
isochronic hypersurface a: t = const (with Ta-{r = 0) = T^h), that guarantees the correctness 
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n 



of the matching (see (l4l. Ilq for details). The analysis provided in Ref. [15|] shows a fairly 
small difference of the one- and two- particle spectra obtained in hHKM and in the case 
of the direct matching of hydrodynamics and UrQMD cascade at the chemical freeze-out 
hypersurface. Thus, in this article we utilize just the latter simplified "hybrid" variant for 
the afterburner stage. 

Let us try to apply the above hydrokinetic picture to the LHC p + p collisions at -y/s = 7 
TeV aiming to get the minimal interferometry radii/volume at the given multiplicity bin. 
As it is known [IJ] the maximal average velocity gradient, and so the minimal homogeneity 
lengths can be reached for a Gaussian-like initial energy density profile. For the same aim we 
use the minimal transverse scale in ultra-high energy p+p collision, close to the size of gluon 



spots 20| in a proton moving with a speed v ^ c. In detail, the initial boost-invariant tube 



^ According to Particle Data Group compilation 



18|. 



ioT p + p collisions has a Gaussian energy density distribution in the transverse plane ej(r 



with width (rnis) R = 0.3 fm 20| and, following Ref. 15|, we attribute it to a initial proper 
time To = 0.1 fm/c. At this time there is no initial transverse collective flow. The maximal 
initial energy density is defined by all charged particle multiplicity bin. The maximum 
initial energy density, ej(r = 0), is determined in HKM, for selected experimental bins in 
multiplicity, by fitting of the mean charged particle multiplicity in that bin. 

The correlation function for bosons in the UrQMD event generator is calculated according 
to: 

Zl ^A(g - Pi + Pi)(l + COs{pj - Pi){Xj - Xi)) 

where Sa{x) = 1 if |x| < Ap/2 and otherwise, with Ap being the bin size in histograms. 
The output UrQMD 3D correlation histograms in the LCMS for different relative momenta 
q = P]^ — P2 are fitted with Gaussians at each kx = ^^^^ bin 

C(q) = 1 + A • exp{-Rl^,ql^^ - i?Le4de - KngQlng)- (2) 

The interferometry radii Rout{kT), Rsideikx), RiongikT) and the suppression parameter A are 
extracted from this fit. 

In Fig{T]we demonstrate the results from hydrokinetic model for the pion interferometric 
radii, comparing them with the ones measured by the ALICE Collaboration at the LHC 21 1 
mp + p collisions at the energy ^/s = 7 TeV. As one can see there is a significant systematic 
overestimate of the predicted interferometry volume Vint = RoutRsideRiong iia p + p collision 
even at the minimal homogeneity lengths possible for the given multiplicity classes. This is 
consistent with the results of the first paper devoted to the same topic "Pion interferometry 



testing the validity of hydrodjTiamical models" 22| . In what follows we shall try to improve 
the results of the semi-classical HKM event generator by means of the quantum corrections 
to them [12 1 . 



III. THE QUANTUM CORRECTIONS TO THE HYDROKINETIC RESULTS 



In 12| it is shown that, for small systems formed in particle collisions (e.g. pp, e^e~) where 
the observed interferometry radii are about 1-2 fm or smaller, the uncertainty principle 
doesn't allow one to distinguish completely between individual emission points. Also the 



phases of closely emitted wave packets are mutually coherent. All that is taken into account 
in the formalism of partially coherent phases in the amplitudes of closely spaced individual 
emitters. The measure of distinguishability and partial coherence is then the overlap integral 
of the two emitted wave packets. In thermal systems the role of the corresponding coherence 
length is played by the thermal de Broglie wavelength that defines also the size of a single 
emitter. The Monte-Carlo method ([1]) cannot account for such effects since it deals with 
classical particles and point-like emitters (points of the particle's last collision). The classical 
probabilities are summarized according to the event generator method ([1]), while in the 
quantum approach a superposition of partially coherent amplitudes, associated with different 



possible emission points, serves as the input for further calculations [12|. Such an approach 
leads to a reduction of the interferometry radii as compared to Eq. ([T]). In addition, the 
ascription of the factor 1 + cos{xi — X2){pi — P2) to the weight of the pion pair in ([1]) 
is not correct for very closely located points Xi and X2 because there is no Bose-Einstein 



enhancement if the two identical bosons are emitted from the same point [12|, |23|. The 
effect is small for large systems with large number of independent emitters. For small 
systems, however, it can be significant and one has to exclude unphysical contributions 



("double counting" [12|) in the two-particle emission amplitude. Such corrections lead to 
a suppression of the Bose-Einstein correlations that is manifested in a reduction of the 
observed correlation function intercept compared with one in the standard method ([1]). 



The results of Ref. 12|] are presented in the non-relativistic approximation related to the 
rest frame of the source moving with four- velocity m'^. In the hydrodynamic/hydrokinetic 
approach the role of such a source at a given pair's half-momentum bin near some value p 
is played by the fluid element or piece of the matter with the size equal to the homogeneity 







length X{p) [7|. These lengths are extracted from the HKM simulations, namely, from the 
interferometry radii defined by the Gaussian fits to the correlation functions obtained in 
HKM. All the pairs in procedure (1) are considered in the longitudinally co-moving system 
(LCMS) that in the boost-invariant approximation automatically selects the longitudinal 
rest frame of the source and longitudinal homogeneity length in this frame (it is Lorentz- 
dilated as compared to one in the global system |5|). The femtoscopy analysis is typically 
related to a fixed pt bin and so one needs also to determine the transverse source size in 
the transverse rest frame. The corresponding Lorentz transformations do not change the 
szde- homogeneity length; as for the oit^direction we proceed in the way proposed in Ref. 5[. 
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FIG. 1. The pion interferometry volume dependency on the charged particles multiplicity at 



Pt = 0.2 — 0.3. The LHC ALICE data [21[ are compared with pure HKM results (blue solid line) 



and with the quantum corrected ones (red lines). 



Remaining within the Gaussian approximation, reahzed for expanding inhomogeneous 
systems in the saddle point method [7|, ISj, let us fix some p = (pr,0) in the basic LCMS 
reference system (basic-KS) and select the transversely moving reference systems (marked 
by the sign tilde) where the emission density distribution, related to the space-time center 
of this local source xoi{p), io^p), can be well approximated as the following 



p(x,t)oce-^«^»'/2^?(P)-*"/2^'(^) 



(3) 



Then in tilde-RS the correlation function has the form ([2]) where Rf^ng = -^Lng' ^tide — -^side' 



p2 _, p2 



X 



out + ^k^T^ with T defining the duration of emission in this tilde-KS. The 
pair's half-momentum p corresponds to the concrete experimental bin taken in the basic-KS, 
the difference of the particle momenta q components in selected tilde-RS are (jout, (jside = Qside, 
qiong = qiong- Therefore, only Qout and correspondingly Rout, including Xout and T, are really 
transformed in ([2]) at the Lorentz boosts along the transverse momentum of the pair. The 
correlation function C{p,q) is the Lorentz invariant, therefore Rlutip)'iout — i^^- 

To relate the interferometry radius in the rest frame of the source (marked by the asterisk) 



to the one in basic-KS one should express both values through the radius in tilde-KS using 

the invariance property similar as it is done in J5|. Then one can get 

cosh yT 



^outVP) ~ ^O'^tyP) „„„uf„. „ N' ^side ~ ^side, ^long ~ ^long (4j 



^* _^ COSht/T Plutrj.* _rj. Sluh |/T 



^ cosh{yT - Vt) ' P*o cosh{yT - Vt) 

Here rjT is a rapidity of the source in transverse direction, y^ = {yir + 1/2t)/2 is half-sum of 

transverse rapidities of the particles forming the pair. Then one can represent the correlation 

function again in the form ([21) where all the variables are related already to the rest frame of 

the source and the HBT radii in this rest frame are expressed through the radii in basic-KS 

according to (jl]). Note that yr = Ht — Vt, and if the rapidity of the pair is equal to the 

rapidity of the source, y^ = 0, then in this particular case the radius in the rest frame is 

Lorentz dilated by the factor 7. Generally, the reference system where the pair's momentum 

is zero does not coincide with the rest frame of the source that emits the pair. Therefore, the 

direct application of these formulas is not an easy task for the rather complicated emission 

structure in a hypothetical hydrodynamic/hydrokinetic model oi p + p collisions. In Figj2] 

one can see the structure of the chemical freeze-out hypersurface with the values of collective 

velocity labeled for high multiplicity p + p events in comparison with the ones for central 

Pb + Pb collisions at LHCcl. Of course, the details of the transformation will be different for 

a string-based event generator, therefore we present the analysis for the radii transformation 

just in the two limiting cases -R*^^ = Rout and R*^^ = 'jRout (7 = coshy^). 

We provide the quantum corrections at each pT-hin in the rest frame of the corresponding 

source using Eq. (jl]) and then come back again to the basic-KS. To preserve the previous 

notations one can suppose that the source rest frame coincides with tilde-KS. In what follows 

the tilde and asterisk marks are omitted and all values are related to the source rest frame. 

To account that due to the uncertainty principle the emitters (strictly speaking emitted 

wave packets) have finite sizes ((Ax)^) ~ k~^ {k is the momentum variance of the particle 

radiation) when defining the lengths of coherence, one should at first consider the amplitude 

of the radiation processes and only then make statistical averaging over phases of the wave 



packets using the overlap integral as the coherence measure 12 1. 



^ Note, that the initial maximal energy densities are close in both these processes and the peculiarities of 
the freeze-out hypersurface and velocity profile in p + p case are caused by the very large gradients of 
initial density because of the small initial transverse size. 
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FIG. 2. The chemical freeze-out hypersurface in HKM in the transverse plane for y/s = 7 TeV p+p 



collisions at 



drj 



17.9 (red line) in comparison with the analogous result scaled in both r, rx 



coordinates by the factor 1/3 for y^ = 2.76 TeV Pb + Pb collisions (blue line). The maximal initial 
energy densities are close in both cases. The maximal velocities are marked in the corresponding 
points on the curves. 



Following 12| we present the quantum state ipxXPi t) corresponding to a boson with mass 
m emitted at the time ti from the point Xj as a wave packet with momentum variance k 
which then propagates freely: 



^.,(p,t) = e^^^'-^^V^(^')/(p) 



(5) 



where ip{xi) is some phase and / defines the primary momentum spectrum /(p) that we 
take in the Gaussian form, 

1 



/(P) = /IP) 



p 

-e 2p 



(6) 



(27rA;2)3/2 

with the variance k"^ = mT. The effective temperature of particle emission in the local rest 
frames in HKM, T, is close to the chemical freeze-out temperature Teh- 

The amplitude of the single-particle radiation from some 4-volume can be written at 



very large times too as a superposition of the wave functions ipxXp) with coefficients p(xj) = 
\/p{xi) that leads in the case of completely random phases to the emitter distribution ([3]) 
in the local rest frames of the sources: 

A{p,t) = c d^XiiJxXP^t)p{xi)^ (7) 

where c is the normalization constant. 

The single- and two-particle spectra, averaged over the ensemble of emission events with 
partially correlated phases (p{x) are 

W{p} = c' /"d^a;ciVe*P("-"')p(x)p(x')/(p)(e*(^(")-^("')) 



W{pi,p2) = c^ /"rf^Xirf^X2(^VirfV2e^(Pi^i+P2^2-^^^i-P2^2) ■ 



The phase averages are associated with corresponding overlap integrals 12 1 



dhtPx{t,r)ijl,{t,T) 



(9) 



2 

where ^x^{t,r) = I3/2 J f{p)e~^^^^~^^^e~^^^^*^~^^d^p are the wave functions of single 
bosonic states in coordinate representation. 

Then the correlation function C(p, q) can be expressed through the homogeneity lengths 
in the local rest frame Rl = ^ongip)^ ^s = KideiP)^ ^o = Kutip) ^^at are expressed 
through the HBT radii obtained from the Gaussian fit ([2]) of the HKM correlation functions 
and transformation law (HI) as described above. 



W{pi,p2) 



W{pi)W{p2) 

_„2 n2 _"k)£o__„2 r2 _i*I)£s__„2 r2 '"'Q-^L (q-P)"'^" 4fc^T^ 

1 + e '° "^^^^ «s«si^3^ iL^Lj^^^^ ^2 1+4*=^^^ _Crf(p,q), (11) 



fe^fl? (q.p)2T2 ^^^2^2 



where k^ = k^/{l + ak'^T'^/m'^), parameter a{k'^R^) is defined from the model numerically 



(it is the order of unity for i? ~ 1 f m and tends to zero for the large sources - see 1^ for 
details), and the subtracted term 

24^0^0 (l+^fep^o) 24fcg4(l+8fcg4) 292fcgfl4(i+8fcgfl|) 

Cd(P,q)=e ('+*'0«o)(l+«'=0«O+«'=0«o) (l+4fcgfl|)(l+8feg4+8fe4i?.|) (l+4fc2fl2^) (^i+8fc2fl|+g,4fl4 ) _ 

2fe2T4{p-q)2 (l+Sfc^T^) 

10 



.2Z.2 ,2^2^ _ iko (l + ik'T^) {l + ^klRD 



^ °^ ' ' ^ kil + 8A;2T2 + 8k*T^) (1 + Sk^Rl + Sk^R^o) 

{1 + SP^Rl + Sk^R'^s) {1 + Sk^Rl + Sk^Ri) J ^ ' 

corresponds to the elimination of the double counting. 

Now we can see that the apparent interferometry radii extracted from the Gaussian fits 
to the correlation function fITT]) are reduced as compared to those obtained in the standard 
approach. 

Particularly, if we neglect the double counting effects, truncate the subtracted term 
Crf(p, q) in (fTTj) . and fit the correlation function with the Gaussian ([2]), we obtain the 
femtoscopic radii Rouu Rside, Riong related to the standard ones Rout,st, Rside,st, Riong,st as 
follows 

^out ( d2 ^k^Rp rr2„.2 ^K T 



/?2 

-^ ''out.st 



(2 O'^O m2 2 \ / / T~i2 m'2 2 \ 



R2 /lt.2D2 

-"-side ^'^0-"'S 



RLg _ AklRl 
Rln,.st 1 + "^klRl 



(13) 



where Vout = Plut/Po ^ ^ according to the non-relativistic approximation. For large source 
sizes, e.g. when the homogeneity lengths correspond to A + A-collisions, k^R^ ^ 1, k'^T'^ ^ 
1, all these ratios tend to unity. 

The mean emission duration is supposed to be proportional to the average system size, 
T = a{Ro + Rs + Rhj/'i that leads to a quadratic equation expressing Rq (and T) through 
Ri^st- The latter are connected with ones taken in the hasic-H^ according to transformation 
laws dl]). The value a is a free model parameter. Then we put these extracted values into 
the expression (ITT]) for the correlation function and perform its fitting with the Gaussian (|2]). 
This gives us finally the interferometry radii Rout-, Rside and Riong in view of the uncertainty 
principle. The radii are presented then in the basic-KS using the transformations inverse 
to (H. 



The correlation function is the ratio of the two- and one-particle spectra. It is found 12 1 
that quantum corrections to this ratio are not so sensitive to different forms of the wave 
packets as the spectra itself. In particular, the effective temperature of the corrected trans- 
verse spectra depends on whether the parameter of mean particle momentum is included 

11 



or not into the wave packet formalism. If yes, the corrected effective temperature for small 
sources i? ~ 1 fm is equal or even higher than that of individual emitters, T = /c^/ttt,, while 



for the wave packets in the form (^ it is lower [12|. Besides of this, in the nonrelativistic 
approximation one can describe only very soft part of the spectra. That is why we focus 
in the article on the corrections to the Bose-Einstein correlation functions where in the rest 
frame of the source the total and relative momenta of the boson pairs are fairly small. 



IV. THE RESULTS FOR p + p AND p + Pb COLLISIONS, AND DISCUSSION 

The initial conditions for HKM are described in Section 2. The HKM event generator 
provides us with the interferometry radii in basic-KS. To find the corresponding homogeneity 
lengths in the rest frame of the source according to (j4]) we use, as discussed in Section 3, 
the two limiting cases: the transverse boost to the rest frame of the pair from the basic 
LCMS system, or no transformation at all. For the former case it is defined by the px bin 
and for p^ = 0.2 — 0.3 GeV coshy^- = ^ = 2.05. The parameter a connecting T with Ri 
increases linearly with multiplicity from 0.67 to 1.0, primary momentum spectrum dispersion 
k = 0.16 GeV/c (T=0.18 GeV), pair mean transverse momentum in the source rest frame 
p^ = 0.1 GeV/c. The a parameter is set to linearly decrease with multiplicity from 1.3 to 0.6 



which is close to the theoretical results [12|. As for the 7 = 1 case, the parameter a decreases 
linearly with multiplicity from 0.65 to 0.55, A; = 0.18 GeV/c and p^ = 0.15 GeV/c. The a 
parameter decreases linearly with multiplicity from 1.2 to 0.78. In Fig. [T] along with the 
experimental and pure HKM results we present the multiplicity dependence of the quantum 
corrected interferometry volume at pr = 0.25 GeV/c. The solid line represents the corrected 
values calculated under the assumption that the Rout interferometry radii, observed in basic- 
KS, are Lorentz-contracted by a factor 7 = 2.05 for the chosen pt = 0.25 GeV/c value as 
compared to ones in the source rest system. The dashed line demonstrates the no-contraction 
case when 7 = 1. As one can see the accounting for the uncertainty principle allows one to 
describe the overall multiplicity dependence of the interferometric radii. Figure [3] represents 
the dependence on multiplicity of individual radius parameters. 

The suppression of the Bose-Einstein correlations for small sources with closely located 
emitters takes place even without specific coherence mechanism and resonance contributions. 
To see this effect the double counting in the correlation function should be eliminated as 
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FIG. 3. The pion interferometry radii dependency on charged particles multipUcity 
at pt = 0.2 — 0.3 GeV. The designations are the same as at the Fig. [TJ 

Eq. ( TTT]) demonstrates. Then the additional suppression parameter X^oh < 1 in the Gaussian 
fit appears and the relevant parameter in ([2]) becomes A = Aco/iAhkm- The result of our 
calculations gives Xcoh = 0.9 - 0.95 for not very small multiplicities. 

In addition to the correlation analysis oip+p collisions, let us make the simplest estimates 
and try to predict the HBT radii for p + Ph collisions at the LHC energy ^/s = 5.02 GeV. 
We ignore the possible asymmetry of the hydrodynamic tube in the longitudinal direction 
and present our prediction within hHKM for centrality c = — 20 % with dNch/dr] = 35. 
The results are calculated for the two initial radii with rms equal to 0.9 fm and 1.5 fm and 
also for the two initial times: r = 0.1 fm/c and 0.25 fm/c. It turns out that the latter factor 
is not essential if we keep fixed final multiplicity: only the longitudinal radii are 3-4% higher 
at r = 0.25 fm/c than at 0.1 fm/c. The transverse radii practically coincide. Therefore, we 
finally demonstrate only the case r = 0.1 fm/c. The initial transverse sizes of the system, 
created in the p + Pb collision, are taken from the Ref. 2J]: "In the conventional wounded 
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nucleoli model it is assumed that the sources are located in the transverse plane in the 
centers of the participating nucleons. This amounts to rather large initial transverse sizes 
in the p — Pb system, R = 1.5 fm. Locating the source in the center-of-mass of the A^A^ 
system is also admissible, which leads to a more compact initial distribution, R = 0.9 fm". 
The results for the interferometry volume are presented at the Fig. |H The model parameter 
set is extrapolated from the described above and consistent with that for the p + p system, 
with 7 = 1 to the case of larger sizes typical for the p + Pb collisions. At that the k and p^ 
values are left the same as for the p + p case, whereas a and a are chosen to be smaller. For 
the R = 0.9 fm initial transverse size a = 0.5, a = 0.43 and for R = 1.5 fm we put a = 0.45, 
a = 0.37. 

Considering the multiplicity dependence of HBT radii in p + p and p + Pb collisions we 
cannot bypass the scaling hypothesis issue 25|, that suggests a universal linear dependence 
of the HBT volume on the particle multiplicity. It means that the observed interferometry 
volume depends roughly only on the multiplicity of particles produced in collision, but not 



on the geometrical characteristics o 



the collision process. At the same time, as it was found 



in the theoretical analysis in Ref. [26[, the interferometry volume should depend not only 
on the multiplicity, but also on the initial size of colliding systems. The figure |4] shows the 
dependency Vint{{dNch/dri)) for the case oi p + p collisions at the LHC, y/s = 7 TeV, and for 
the most central (only!) collisions of nuclei having similar sizes, Pb + Pb and Au + Au, at the 
SPS, RHIC and LHC. We have also added on the plot our prediction for the interferometry 
volume of pP6-system, that has an initial size larger than that for the pp-system. As one can 
see, the different groups of points corresponding to p + p, p + Pb and A + A events cannot 



be fitted by the same straight line. This apparently confirms the result obtained in 26|] that 
the interferometry volume is a function of both variables: the multiplicity and the initial 
size of colliding system. The latter depends on the atomic number A of colliding objects 
and the collision centrality c. 

At the figures [5] and [6] for the two multiplicity classes {dNch/drj) = 9.2 and {dNch/drj) = 
17.9 we present the three curves for interferometry radii as a function of px'- the experi- 
mental one, the one taken just from the HKM simulations and the other one obtained after 
application of the quantum corrections. The basic parameters used correspond to the case 
7 = 1 (see above). The a parameter values linearly increase with px from 1.0 to 1.1 for the 
{dNch/drj) = 9.2 case and from 0.85 to 1 for {dNch/dr]) = 17.9 in such a way that for p^-bin 
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FIG. 4. The interferometry volume dependency on charged particles multiplicity. The curve 
fragments in the middle correspond to the HKM prediction for p + Pb collision at the LHC energy 
^/s = 5.02 GeV. The upper one is related to the initial transverse system size R = 1.5 fm and for 
the two lower ones R = 0.9 fm. The curves at the left and at the right represent the HKM results 



for p+p and j4 + A central collisions respectively. 
RHIC and LHC, taken from papers [21[, [23] - 



compared to the experimental data at AGS, SPS, 



34l |. The pp volumes are calculated as a product 



RoutRsideRiong of respective experimental radii. The blue lines correspond to pure HKM results, 
whereas the quantum corrections to them are presented by the red lines. 

(0.2,0.3) it has the same values as in previous Ri{{dNch/dri)) calculations. As one can see, 
similarly to the multiplicity behaviour, the quantum corrected p^-dependency of the radii 
also gets closer to the experimental values, but for large pt the corrections are insufficient 
to fully describe the observable femtoscopy scales behavior. This fact may indicate that 
sources of particles with large px cannot be described in hydrodjTiamical approximation. 
Note that just for such large pr the non-trivial base-line corrections, already provided in 
presented experimental data, are very essential. 

Finally we demonstrate the prediction of the pT-dependence of the radii oi p + Pb colli- 
sions. It is presented in FigJT] together with the corrections due to the uncertainty principle. 
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FIG. 6. Dependence of interferometry radii on p^, {dNch/drj) = 17.9 
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FIG. 7. The HKM prediction for the dependence of p + Pb interferometry radii on pj- at the LHC 
energy ^ = 5.02 GeV, [dN^h/dr]) = 35. 



One can see that corrections are smaller for the systems with larger homogeneity length 
and not very essential for p + Ph and probably ioi A + A peripheral collisions. As for a 
homogeneity lengths formation, in a very recent paper 35|] it is found that, in the absence 
of hydrodynamic flow, HBT radii should be similar in pp and pPb collisions. 



V. CONCLUSIONS 



One can conclude that quantum corrections to pion interferometry radii m p + p col- 
lisions at the LHC can significantly improve the (semi-classical) event generator results 
that typically give an overestimate of the experimental interferometry radii and volumes. 
The corrections account for the basic (partial) indistinguishability and mutual coherence of 
closely located emitters because of the uncertainty principle [12|. The additional suppression 
of the Bose-Einstein correlation function also appears. The effects become important for 
small sources, 1-2 fm or smaller. Such systems cannot be completely random and so require 
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a modification of tlie standard tlieoretical approacli for correlation analysis. The predicted 
interferometric radii for p + Pb collisions need some small corrections only for its minimal 
values corresponding to the initial tiansveise size oi pPb system 0.9 fm. 

More sophisticated result of this study is a good applicability of the hydrodynam- 
ics/hydrokinetics with the quantum corrections for description of HBT radii not only in 
y4 + A collisions but also, at least for large multiplicities, in p + p events. These radii are 
well reproduced for not too large pt- Whether it means the validity of the hydrodynamic 
mechanism for the bulk matter production in the LHC p + p collisions is still an open 
question. It is also related to the problem of early thermalization in the processes of heavy 
ion collisions; the nature of this phenomenon is still a fundamental theoretical issue. 
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